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ABSTRACT 

We present Very Large Array images of a "Double-Double Radio Galaxy", a class 
of objects in which two pairs of lobes are aligned either side of the nucleus. In this 
object, B0925+420, we discover a third pair of lobes, close to the core and again 
in alignment with the other lobes. This first-known "Triple-Double" object strongly 
increases the likelihood that these lobes represent mutiple episodes of jet activity, as 
opposed to knots in an underlying jet. We model the lobes in terms of their dynamical 
evolution. We find that the inner pair of lobes is consistent with the outer pair having 
been displaced buoyantly by the ambient medium. The middle pair of lobes is more 
problematic - to the extent where an alternative model interpreting the middle and 
inner "lobes" as additional bow shocks within the outer lobes may be more appropriate 
- and we discuss the implications of this on our understanding of the density of the 
ambient medium. 
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1 INTRODUCTION 

The behaviour of FanarofT-Riley II (FRII; Fanaroff & Riley 
1974) radio galaxies is typified by the ejection of highly col- 
limated, laminar jets from a central core. These jets travel 
supersonically away from the core and are decelerated by 
the ram- pressure of the intergalactic medium (IGM). This 
deceleration results in a bow shock being driven through 
the IGM, inflating large "cocoons" of low-density material, 
part of which we observe as radio-lobes. Finally the energy 
is lost in bright hotspots of shocked emission. Although the 
intrinsic evolution of these powerful radio sources is now 
largely understood (e.g. Kaiser & Alexander 1997, hereafter 
KA; Blunder!, Rawlings & Willott 1999), there are still sev- 
eral unanswered fundamental questions relating to the duty- 
cycle of the active galactic nucleus (AGN) and the contents 
of the lobes. 

"Double-Double Radio Galaxies" (DDRG) are FRIIs 
with a second pair of lobes; all four lobes appear to em- 
anate from a common nucleus and are strongly coaligned to 
within 10° of each other (although see also Saikia, Konar & 
Kulkarni (2006), who identify other DDRG whose lobes are 
misaligned). When viewed in isolation each of these lobes 
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has a comparable morphology with that of a "normal" FRII 
galaxy. A sample of seven such objects has been identified, 
studied and modelled by Schoenmakers et al. (2000a, b), 
Kaiser, Schoenmakers & Rottgering (2000) and Konar et 
al. (2006); these seven objects are B0925+420, B1240+389, 
B1450+333, B1834+620, 3C219, 3C445 and 4C26.35 (al- 
though in retrospect the inner structures of 3C219 are very 
narrow and may instead be "knots" in the jet rather than 
new lobes). Saikia et al. (2006) studied a larger sample of 
DDRG, bringing the total to twelve. The additional pair of 
radio lobes has been interpreted as evidence for a second 
phase of jet activity, the inner pair of lobes forming within 
the low-density cocoon of the original jet. With this unusual 
opportunity to obtain information about multiple episodes 
of activity, this subset of the FRII class has the potential to 
give us rare insight into the conditions of the surrounding 
medium before and after the launching of the jet, as well as 
duty-cycle timescales. 

There are still uncertainties regarding certain aspects of 
these objects. Are the inner structures truly separate radio 
sources with no physical connection to the outer lobes? Are 
the ages and ambient densities of the inner lobes - which 
may be the density of the cocoon inflated by the original jet 
- atypically low, i.e. much lower than these values are for 
objects of similar size as the inner radio lobe? Are there any 
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Table 1. Angular resolutions, position angles (PA) and la noise levels in each of the new VLA images. 



Frequency 
(GHz) 


A-array Resolution (PA) 
(") (°) 


A-array Sensitivity 
mjy beam -1 


B-array Resolution (PA) 
(") (°) 


B-array Sensitivity 
mjy beam -1 


1.4* 






10 x 10(0) 


3.21 X 10 _1 


1.4 


1.7 x 1.5(5.1) 


3.15 X 10~ 2 


6.4 x 5.0(68.2) 


1.70 X 10" 1 


4.9 


0.6 x 0.5 (-28.7) 


4.28 x 10~ 2 


1.7 x 1.5 (-68.8) 


2.72 x 10~ 2 


8.4 


0.3 x 0.2 (-23.5) 


3.20 x 10~ 2 


0.9 x 0.8 (-68.5) 


2.22 x 10~ 2 



*The new B-array VLA data, combined with the FIRST data and convolved with a 10" beam as displayed on the left-hand side of Fig. 1. 
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Figure 1. VLA images of B0925+420 showing the three pairs of lobes. All contours are at -3, 3, 6, 12, 24, 48x<r, with the exception of 
the left-hand plot which has an additional 2a contour. LHS: the 2001 VLA 1.4-GHz B-array data, combined with additional 1.4-GHz 
B-array data from the 1994 FIRST survey (Becker, White & Helfand 1995) and finally reconvolved with a 10" beam in order to detect 
the weak outer lobes at a 3a level (a = 3.21 X 10 _1 mjy beam -1 ). Centre: 2001 VLA 4.9-GHz B-array image, clearly resolving the middle 
and inner lobes but resolving out all emission from the outer lobes (a = 2.72 X 10~ 2 mjy beam -1 ; angular resolution = 1.7 X 1.5"). RHS: 
2001 VLA 8.4-GHz B-array image showing a close-up of the inner lobes (a = 2.22 X 10~ 2 mjybcam -1 ; angular resolution = 0.9 X 0.8"). 



remnants of jets and hotspots in the outer lobes? Is the den- 
sity much lower than in the "normal" galactic environment, 
which would strongly suggest that the inner structure is in- 
deed embedded in the "hole" excavated by the outer lobes? 
Or is the density higher which would imply very efficient re- 
placement of the inner lobes by the heavier external medium. 
Either result would be extremely important because of the 
constraints placed on the physics. 

In this paper we study one such DDRG, B0925+420 
(=J0929+4146), which provides particularly strong evidence 
for multiple episodes of activity. We present the observations 
in Section 2, analysis of the images in Section 3 and details 
of the polarisation properties in Section 4. In Section 5 we 
present the results of analytical modelling of each of the ra- 
dio lobes. We make comparisons with other radio transients 
in Section 6 and finally draw our conclusions in Section 7. 

2 OBSERVATIONS 

Radio observations of B0925+420 were obtained using the 
VLA on 2000 October 20-21, when the array was in the 



highest resolution A-configuration, and 2001 April 2, when 
the array was in the B-configuration. Data were obtained 
at three frequencies - 1.4-GHz (L-band), 4.9-GHz (C-band) 
and 8.4-GHz (X-band) - in each array; IFs at two adjacent 
frequencies were used in each continuum observation. In or- 
der to reduce potential bandwidth-smearing in the more ex- 
tended images, all A-array observations and the 1.4-GHz B- 
array observations were obtained in the PA/PB spectral line 
observing mode. This gave us 8 channels of 3.125 MHz each, 
resulting in a 25-MHz bandwidth instead of the more usual, 
continuum bandwidth totalling 50 MHz. Full polarisation in- 
formation was also recorded. 3C 286 (J1331+305) was used 
as the flux and polarisation calibrator while JVAS 0920+446 
was used for phase calibration and to ensure parallactic an- 
gle coverage. 

Additional images were obtained from the Faint Im- 
ages of the Radio Sky at 20-cm (FIRST; Becker et al. 1995), 
NRAO VLA Sky Survey (NVSS; Condon et al. 1998) and 
Westerbork Northern Sky Survey (WENSS; Rengelink et al. 
1997) surveys at 1.4-GHz, 1.4-GHz and 0.33-GHz respec- 
tively. 
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Figure 2. Spectral index maps created from the 2001 VLA data; we used only the B-array images so as to ensure that no emission was 
resolved out at the higher frequencies. In each case, the higher-resolution image was reconvolved with the beam of the lower-resolution 
image. For comparison with the spatial scale of the original high-resolution image, we overplot each map with the 3cr Stokes I flux density 
as measured in the unconvolved image. LHS: 4.9-GHz/8.4-GHz with angular resolution = 1.7 X 1.5". Inset plots show the northern and 
southern middle lobes more clearly and we note that the regions of unexpectedly steep spectral index correspond to those regions outside 
the 8.4-GHz contour. RHS: 1.4-GHz/4.9-GHz with angular resolution = 6.4 X 5.0". The core appears to have a relatively fiat spectrum 
and the lobes a steeper spectrum with a in the range —1 to —2; the spectrum of the middle lobes appears steeper than that of the inner 
lobes. 



The 2000/2001 VLA data were reduced using standard 
flagging, calibration, cleaning and imaging routines within 
MPS. The flux densities of 3C 286 were obtained using the 
formulae of Baars et al. (1977) but with the revised coeffi- 
cients of Rick Perley, as is the default option in the SETJY 
routine. Additional analysis routines within MPS were used 
to convolve the images with different beam-sizes, combine 
the new 1.4-GHz data with the FIRST survey image, mea- 
sure lobe parameters and create spectral index maps. We 
list the angular resolutions and sensitivities of the images 
at all three frequencies and in both array-configurations in 
Table 1. 



3 RESULTS - IMAGES 

B0925+420 is detected in the 2000/2001 VLA data at each 
of the six frequency /resolution combinations. While much 
of the extended emission is resolved out in the higher- 
resolution A-array images, the remaining images clearly 
show a series of four well-aligned FRII-like lobes. Closer in- 
spection reveals that the outer lobes detected in the FIRST 
image (Schoenmakers et al. 2000a) are resolved out in the 
more recent images and that the core has been resolved at 
high frequency into an additional, third, pair of lobes. 



We have obtained the 1.4-GHz FIRST image and com- 
bined the data with the more recent 1.4-GHz data. Finally, 
by convolving this resultant image with a larger beam (10") 
we are able to detect the outer lobes at a 3<r level. This 
reconvolved image is shown on the left-hand side of Fig. 1. 
We note that the morphology of the outer lobes is consistent 
with that shown in Schoenmakers et al. (2000a); it is sim- 
ply the choice of contour levels which leads to this apparent 
inconsistency. The second pair of lobes is resolved clearly 
while the core remains unresolved. The source to the left 
of the image is unrelated to B0925+420. Comparison with 
the NVSS image suggests that no flux is resolved out of the 
middle lobes in the new VLA B-array images at 1.4-GHz. 

The middle of Fig. 1 is the 4.9-GHz B-array image. The 
outer lobes are resolved out but the middle pair of lobes 
are shown clearly. The core is now resolved into a smaller 
pair of lobes. Finally, the right-hand side of Fig. 1 shows the 
8.4-GHz B-array image; even the inner lobes are resolved 
out and/or undetected in the A-array image. A bright core 
is still detected, suggesting the possibility that either these 
lobes are still being fed by jets or that the core may yet be 
further resolved into more components. 

Following the arguments outlined above for the "double- 
double radio galaxies" providing the best-known evidence 
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Table 2. Dimensions and flux densities (S v ) of the lobes obtained from the images. Typically the lower-resolution images were used so 
as to reduce the likelihood of resolving out any flux. Errors are to within up to ~ 5" and ~ 1 mjy. Lobe-lengths have been measured 
from the leading edge to the centre of the core; widths have been measured at their widest points in the non-convolved images. The last 
two columns list the spectral indices in the 1.4-GHz / 4.9-GHz and 4.9-GHz / 8.4-GHz bands. These values are calculated from the flux 
densities and arc broadly comparable with the spectral index maps shown in Fig. 2; we note that (i) the calculated spectral index for the 
southern middle lobe is artificially steep due to some flux being resolved out at 8.4-GHz, and (ii) the calculated spectral index of —1.22 
for the core is artificially steep due to contamination by emission from the inner lobes at low resolution. 



Lobe 


Length 

(arcsec) 


Width 

(arcsec) 


S 0.32 MHz 
(mjy) 


Si. 4 GHz 
(mjy) 


S 4.9GHz 
(mjy) 


5 8.4GHz 
(mjy) 


a 

(1.4-GHz/4.9-GHz) 


a 

(4.9-GHz/8.4-GHz) 


North, outer 


227 


38 


88 


23 










South, outer 


193 


42 


84 


24 










North, middle 


83 


4 




32 


10 


4 


-0.93 


-1.70 


South, middle 


57 


4 




29 


10 


2 


-0.85 


-2.99 


North, inner 


4 


2 






1.4 


1.0 




-0.62 


South, inner 


3 


1 






0.5 


0.3 




-0.95 


Core 








6 


1.3 


1.4 


-1.22 


0.14 



% 



% 
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Figure 3. Fractional linear polarisation (FP) image at 1.4-GHz 
using the 2001 VLA:B-array data only; FP is represented by the 
grcyscale. The vectors indicate the direction of the electric field, 
following correction for RM, and their length indicates the linear 
polarisation (1" = 4.46 X 10" 5 Jy beam" 1 ). The core and inner 
lobes are weakly polarised; the middle lobes are polarised at the 
~ 10 — 20% level, higher in the southern lobe and around the 
lobe edges. The position angle of the vectors indicates that the 
magnetic field in the middle lobes is aligned with the long axis of 
the lobe, although there appears to be some rotation in the inner 
lobes. 



for the duty cycle of jets in radio galaxies, we propose that 
B0925+420 is therefore the first known radio galaxy 
which has been observed to display three episodes 
of jet activity. This "triple-double" radio galaxy provides 
us with the strongest evidence yet that these objects are 



indeed potential laboratories for the study of the duty cy- 
cle of AGN jets, as opposed to a sequence of bright knots 
in an underlying jet structure. We discuss this suggestion 
quantitatively in Section 5. 

Finally we have created spectral index maps from the 
1.4-GHz / 4.9-GHz and 4.9-GHz / 8.4-GHz images; we define 
the spectral index, a, in terms of S v oc v a , where S v is the 
flux density at frequency v. In order to reduce the chances of 
flux being resolved out, only the B array images were used. 
In each case, the higher- frequency image was reconvolved 
with the beam of the lower-frequency image and the resul- 
tant images presented in Fig. 2. The core shows a very flat 
spectrum with a close to zero. The lobes show steeper spec- 
tra of —1 to —2, with the middle lobes appearing to have 
steeper spectra than the inner lobes. The spectral index of 
the southern middle lobe, in particular, appears artificially 
steep. If we assume a value of cn.g-GHz/s,A-GHz = —2, it 
suggests that the expected flux density at 8.4-GHz is ~ 3 — 4 
mjy and so ~ 2 mjy of flux has been resolved out; sim- 
ilarly, assuming a spectral index of -1 suggests that ~ 5 
mjy has been resolved out. These estimates require that 
no flux is resolved out at 4.9-GHz, a reasonable assump- 
tion given the uniformity of the southern middle lobe in the 

Ql.4-GHz/4.9-GHz map. 

While the bright core suggests the possibility that the 
inner lobes are still receiving energy via jets (as also found 
in other objects, e.g. J0041+3224; Saikia et al. 2006), the 
absence of hotspots in the outer and possibly the middle 
lobes suggests that these structures are no longer supplied 
with energy. Such a statement is supported by the steep 
spectral indices presented in Fig. 2 (middle lobes) and the 
table of observed parameters (Table 2), which includes ad- 
ditional data from the Westerbork array and thus provides 
a second frequency at which the outer lobes were detected. 
Theoretical modelling in Section 5 adds further support to 
these observations. We note that we cannot categorically 
rule out the presence of hotspots in the middle lobes - some 
3CR objects have recessed hotspots, and structures of sizes 
comparable with the resolution limit have been detected in 
other objects (e.g. Jeyakumar & Saikia 2000) - but, with no 
firm evidence to the contrary, we assume that they are not 
present in the middle or outer lobes of B0925+420. 
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4 RESULTS - POLARISATION 

The middle and inner lobes of B0925+042 were also de- 
tected in the linear polarisation (LP) images made from the 
2001 VLA data in B array. We do not have polarisation 
calibration available for the FIRST data and so cannot ex- 
tend this analysis to the outer lobes. The LP flux density at 
4.9- and 8.4-GHz was low and/or the structure partially re- 
solved out but the structure is clear at 1.4 GHz and, where 
structure was visible, the fractional polarisation (FP) did 
not appear to vary significantly with frequency. A rotation 
measure (RM) image was created, using all three frequen- 
cies. The lowest possible RM values were in the range 10-15 
radm~ 2 , which is consistent with the Galactic medium be- 
ing the only Faraday screen between us and the source; such 
a conclusion was also reached for B1834+620 (Schoenmak- 
ers et al. 2000b) and any intrinsic contribution to the RM 
in J1453+3308 was deemed to be small (Konar et al. 2006). 
We show the 1.4-GHz FP image in Fig. 3; electric field vec- 
tors are over-plotted and include the correction (28°) due to 
Faraday rotation. 

The middle lobes are clearly polarised to 10-20%, higher 
in the southern lobe and towards the edges; such a FP is 
comparable with that measured for J1453+3308 (Konar et 
al. 2006). The core and/or inner lobes of B0925+480 are 
weakly polarised. The direction of the electric field vec- 
tors appears to change between the two pairs of lobes - 
at the centre, the magnetic field is apparently aligned quasi- 
perpendicularly to the position angle of the radio galaxy and 
direction of jet motion. The magnetic field in the middle 
lobes has become aligned comparably with the radio galaxy. 
We note that this effect in the core/inner lobes is ambiguous 
since they are not fully resolved and we may be looking at 
the combined effects of an unpolarised core and polarised in- 
ner lobes. This structure of the magnetic field in B0925+480 
seems less complex than in J1453+3308, in which the vec- 
tors rotated at the edge of the lobe (Konar et al. 2006); 
however, as discussed in Sections 5, we highly suspect that 
the outer lobes of J 1453+3308 exist in a very different envi- 
ronment from the middle lobes of B0925+480 and so we are 
not comparing like with like. 

The low RM and high FP of B0925+480 suggest that 
there is no evidence for depolarisation in the middle lobes, 
hinting that the environment is not very dense and/or 
clumpy. With their lower levels of LP, the inner lobes could 
still be dense and clumpy. The modelling in the next section 
provides further clues as to the nature of the environment of 
each pair of lobes, confirming the results of the polarisation 
observations. The alignment of the magnetic field with the 
edge of the middle lobes may be indicative of compression 
of the field by the lobe expansion. Surprisingly, given the 
absence of evidence for a current energy supply, this would 
argue for continuing expansion of the middle lobes. 



5 MODELLING THE SOURCE-SIZE AND 
LUMINOSITY 

5.1 Brief description of the model 

We use the model of KA for the dynamical evolution of the 
large-scale structure of extragalactic radio sources with a 
FRIl-type morphology. The model envisages two jets prop- 



agating in opposite directions from a central AGN. The ram 
pressure of the jet material is balanced by the ram pressure 
of the receding external medium in front of the jet flow. This 
interaction between the jet and its environment gives rise to 
a strong shock at the end of the jet and a bow shock in 
the external gas. After passing through the shock, the jet 
material inflates a lobe or cocoon surrounding the jet flow. 

The thermal pressure inside the lobe, p, balances the 
pressure of the jet material and thus confines the jet flow. 
The lobe pressure also drives the sideways expansion of the 
lobe. In the model of KA the lobes are overpressured with 
respect to the external medium and so the sideways expan- 
sion of the lobes is also confined by the ram pressure of the 
receding external medium. In this way, the bow shock in 
front of the jet shock is extended to envelope the entire lobe 
structure. The density distribution in the external medium 
is modelled as a power-law centred on the AGN with 

' = W 

where p and ao are constants and r measures the distance 
from the AGN. The model then predicts the length of the 
lobe in the direction of the jet flow, L, as a function of the 
age of the flow, t, the energy transport rate of the jet or jet 
power, Q, the exponent for the power-law density distribu- 
tion, /3, and the 'density parameter', poa^- Another param- 
eter entering the calculations is the aspect ratio of the lobe, 
R, which we define here as the ratio of L and the width of 
the lobe, W, measured at the widest point along L so that 
R = L/W. 

The model of KA was extended by Kaiser, Dennett- 
Thorpe & Alexander (1997, hereafter KDA) to allow the 
prediction of the radio luminosity density of the lobe, P„, 
as a function of essentially the same parameters. The KDA 
model assumes that the radiating relativistic electrons arc 
injected into the lobe at the jet shock with a power-law 
energy distribution according to n(E) dE oc E~ s dE. The 
subsequent energy losses of the electrons due to the adia- 
batic expansion of the lobe and radiative processes, emission 
of synchrotron radiation and inverse Compton scattering of 
cosmic microwave background photons, are taken into ac- 
count self-consistently. 

The middle and (particularly) the outer lobes of 
B0925+420 show little evidence for the presence of radio 
hotspots which are usually identified as the impact sites of 
the jets and the location for the injection of relativistic elec- 
trons into the lobes. It is therefore likely that the energy 
supply of the jets as well as the injection of relativistic elec- 
trons into these lobes has ceased some time ago. If we assume 
that some disturbance of the jet flow close to the AGN is 
responsible for stopping the jet at a time t s after the onset 
of jet activity, then the injection of relativistic electrons into 
the lobe stops at t B + t t , where tt is the time it takes the last 
jet material to travel from the AGN to the end of the jet. 
Assuming that the jet material travels with a speed close to 
the speed of light inside the jet, we have tt = L/c where L 
is a function of t B + t t . Using the model of K A we can solve 
this implicit relation for tt by iteration. 

When the energy supply of the jet to the lobe shuts 
down at time t B + U, the dynamics of the lobe expansion 
start to change (e.g. Kaiser & Cotter 2002). However, this 
new expansion regime takes a few sound crossing times of the 
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Table 3. Model parameters used in the fitting process. 



Parameter 


Description 


Nature of 


Value if set 






parameter 


ab initio 


u 


Projected length of lobe 


observed 




p 


T ,i i m i Ti(~it;i i"\.r f\r , Yt<zit\r i~it |p»rir» fit - rincrTviTicf froriiinvii'M/ 1/ 


("I icon rn( 1 




Rl 


ilaULL Lr L Cb\j LVJ U 1 ±l_J L/L, 








Source rcdshift 


observed 


0.365 


r 

r c 






4/3 


r B 


Adiabatic index of magnetic field 'fluid' 


set 


4/3 


r x 


Adiabatic index of external medium 


set 


5/3 


a 


Viewing angle of the jet axis to our line of sight 


set 


7I-/2 


k 


Energy density of non-radiating particles in the lobe 


set 





5 


Exponent of initial power-law energy density of relativistic electrons 


set 


2 -► 2.5 


Tmin 


Minimum Lorcntz factor of relativistic electrons at time of injection into the lobe 


set 


1 


7max 


Maximum Lorentz factor at same time 


set 


10 8 





Exponent for the power-law density distribution 


set 


1.5 


CLQ 


Core radius of density distribution of external medium 


set 


2 kpc 


Po 


Density of external medium at a distance ao from the AGN 


fitted 




Q 


Jet power 


fitted 




t 


Time elapsed since jet flow started 


fitted 




t s 


Time at which jet flow is interrupted 


fitted 





Table 4. Input parameters for the model. Angular sizes and flux densities were converted into luminosity densities and lengths from 
the values in Table 2 and using a value of 0.365 for the redshift, z. A WMAP cosmology was assumed, with Ho = 71km s -1 Mpc -1 , 
f2m = 0.27, Ha = 0.73 (Spergel et al. 2003). We note that the luminosity at 1.4-GHz for the outer lobes has been derived from the flux 
density quoted in Schoenmakers et al. (2000a); it appears that some flux was resolved out of the images presented here. 



Lobe 


Length 
(kpc) 
L 


Aspect Ratio 
R 


Freq. 1 
(GHz) 
v\ 


Lum. Density 1 
(10 24 WHz- 1 ) 


Freq. 2 
(GHz) 

^2 


Lum. Density 2 
(10 24 WHz" 1 ) 
Pi/ 2 


North, outer 


1150 


6 


0.32 


25 


1.4 


10 


South, outer 


972 


4.6 


0.32 


24 


1.4 


11 


North, middle 


419 


21 


1.4 


15 


4.8 


4.4 


South, middle 


287 


14 


1.4 


13 


4.8 


4.5 


North, inner 


20.2 


2.6 


4.8 


0.64 


8.5 


0.45 


South, inner 


15.6 


2.8 


4.8 


0.23 


8.5 


0.14 



lobe to establish itself. For the large sizes of the middle and 
outer lobes the time for this adjustment is comparable to 
the total source age. Therefore for simplicity we continue to 
use the model of KA for the expansion of the lobe. However, 
when calculating the luminosity density of the lobe we only 
take into account relativistic electrons injected before the 
time t s + it- The same approach was used in Kaiser et al. 
(2000). 



5.2 Fitting process 

We can use the models of KA and KDA to constrain the 
properties of the jets and the environment of B0925+420 
from the observed properties of its lobes. Table 3 summarises 
the model parameters used in this process. Most of the 
model parameters we do not discuss in detail here, since 
they are described in the papers by KA and KDA. Observed 
parameters are those taken from observations and some of 
these obviously vary from lobe to lobe. Set parameters de- 
scribe the hydrodynamic and other properties of the lobe 
contents and the external medium. They are set to reason- 
able values and are not varied in the fitting process. Fitted 



parameters are the free parameters of the model constrained 
by the observational data. 

We will see below that the very large observed size of 
the outer lobes of B0925+42 imply very old ages for the jet 
flows. If the jets and therefore the lobes are oriented at small 
angles to our line of sight, then the actual lobe length will 
be even larger. Hence we adopt a viewing angle a = 7r/2, 
i.e. the lobes are in the plane of the sky and the actual lobe 
length and aspect ratio equal the observed values. The model 
depends only on the density parameter poof^ rather than on 
ao and po separately. For convenience we set ao = 2 kpc, a 
value appropriate for isolated elliptical galaxies, and j3 = 1.5 
(Fukazawa, Makishima & Ohashi 2004) and only fit the value 
of po. In comparing model results we therefore note that a 
change of po may also imply a change in ao. Finally, the 
slope of the power-law energy distribution of the relativistic 
electrons strongly influences the slope of the observed radio 
spectrum. We adjust the value of 8 in the range 2 to 2.5 as 
discussed in the next sub-section. 

For each individual lobe we use the lobe length and 
aspect ratio as well as the luminosity densities measured at 
two frequencies in the fitting process. Table 4 summarises 
the model inputs derived from our observations. 
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5.3 Model results 



To determine the model parameters po, Q, t and t B , we ran- 
domly choose a large number of combinations of these four 
parameters and calculate the prediction of the model for the 
lobe length L and the luminosity densities P v at two observ- 
ing frequencies. A given parameter combination is deemed 
to be consistent with the observations, if the model results 
are all within 10% of the observed values of L and the two 
values of P v . Each lobe is fitted individually and the results 
for the source age t and the jet power Q are shown in the 
top panel of Fig. 4 where we have set 8 — 2 for all lobes. 

It is re-assuring that for most possible jet powers the 
lobe age reflects the lobe size with the outer lobes the oldest 
and the inner lobes the youngest. We would expect that the 
jets inflating each pair of lobes, outer, middle and inner, 
have the same jet power and age on both sides of the AGN. 
Hence we expect the patches in Fig. 4 (top) for a given 
pair to show at least some overlap. While this is the case 
for the inner and outer lobe pair, the middle pair shows no 
overlap at all. Also, the area in the Q-t plane allowed for 
both lobes of the outer pair is only small. This may be a 
result of a wrong choice for the slope of the initial power- 
law energy spectrum of the relativistic electrons, S, for some 
of the lobes. The radio spectra of the middle lobes are steep 
compared to the spectra of most of the other lobes and may 
imply a steeper initial energy spectrum. In the interests of 
exploring the parameter space, in the middle panel of Fig. 4 
we show the result of changing 8 to 2.5 for the middle lobes. 
We also change 8 to 2.2 for the northern, outer lobe and to 
2.5 for the southern, inner lobe. Both of these also show at 
least slightly steeper radio spectra. This demonstrates the 
effect of changing 8 on the model results. 

The agreement of the model parameters for the two 
outer lobes is now increased, as is that for the inner lobes. 
However, the model is unable to reconcile the parameters of 
the middle lobes, even for values of 8 greater than 2.5. A 
further problem with the middle pair of lobes is illustrated 
in the bottom panel of Fig. 4 which shows the allowed com- 
binations of Q and the density of the external medium, p . 
For a given jet power the model requires the density of the 
gas in the environment of the middle lobes to exceed that of 
the outer lobes. The same result is found for different values 
of 8. The observed source geometry strongly suggests that 
the middle lobes are located in a region where the outer 
lobes have displaced the ambient medium. We argue below 
that at least some of this region has subsequently been buoy- 
antly refilled with ambient gas, but it is very unlikely that 
this refilling would lead to a substantial increase in the gas 
density compared to the undisturbed gas distribution. This 
suggests that the model described here cannot adequately 
fit the observed properties of the middle lobe pair. 

In contrast the inner lobe pair is consistent with a some- 
what reduced gas density compared to the situation for the 
outer lobe pair. Note that asymmetries in the gas distri- 
bution around the host galaxy are very likely on scales of 
several to hundreds of kpc. The absence of any significant 
overlap of the patches for the inner and outer lobe pairs in 
Fig. 4 (bottom) is therefore acceptable. Finally, if the dis- 
ruption of the jet flow to individual lobes is indeed caused 
by processes close to the AGN, then we would expect the 
time is at which the disruption occurs to be the same for 
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Figure 4. Model solutions for all lobes with, in each case, light 
grey patches representing the northern lobe of each pair and the 
dark grey patches representing the southern lobes. Top: allowed 
combinations of jet power and source age for 8 = 2. Middle: al- 
lowed combinations of jet power and source age when we change 
<5 to 2.5 for the middle lobes and the southern, inner lobe, plus 
change <5 to 2.2 for the northern, outer lobe. Bottom: allowed com- 
binations of jet power and the density of the external medium for 
the same values of <5 as in the middle panel. 
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Figure 5. Schematic diagram illustrating the results of the mod- 
elling. The outer lobe, containing the middle lobe, has started 
to rise buoyantly and is being replaced with the heavier ambi- 
ent medium. The darker shading represents the part of the lobe 
which may still be over-pressured with respect to the surrounding 
medium and so may continue to drive a bow-shock. 



both lobes in a given pair. We find that this is the case for 
the inner and outer lobe pairs for the values of 5 discussed 
above. 

We can derive another, independent estimate for the 
power of the currently-active jets. In our model we assume 
that the inner lobes are currently supplied with energy by 
active jets from the AGN. These jets also give rise to the 
radio emission of the unresolved core in our highest resolu- 
tion radio maps. We measure the flux density of the core 
as 6, 1.3 and 1.4 mjy at 1.4, 4.9 and 8.4 GHz respectively. 
The radio spectrum is therefore flat at least at high fre- 
quencies. The greater flux density at the lowest observing 
frequency is probably caused by unresolved flux of the in- 
ner lobes contaminating the core flux. The jet power of the 
currently active jet can be estimated from the flat spectrum 
core flux as 1.1 x 10 37 W (Hording, Fender & Migliari 2006). 
This is consistent with the lower end of the power estimates 
arising from the observed lobes. 

Our findings can be summarised as ranges of model pa- 
rameters consistent with the observed properties of the inner 
and the outer lobe pairs. They are as follows: 

• The power of the jets which inflated the outer lobes was 
in the range 10 38 W to 8 x 10 38 W while the inner lobe pair 
implies a jet power of 10 37 W to 4 x 10 37 W. A similar jet 
power of 1.1 x 10 37 W is estimated from the core emission. 
The jet power appears to have decreased over time. 

• The density at ao in the environment of the outer lobes 
is in the range 8 x 10" 24 kgm" 3 to 8 x 10~ 22 kgm~ 3 . For 
the inner lobes this decreases to 7 x 10 _25 kgm~ 3 up to 
4 x 10~ 23 kgm -3 . 

• The age of the outer lobes lies in the range 25 to 270 Myr 
while that of the inner lobes is 0.4 to 2 Myr. 

• The jet flow supplying the outer lobes with energy is 
predicted to have ceased between 4 and 70 Myr ago. The 
properties of the inner lobes are consistent with a currently 
continuing energy supply through active jets. 



5.4 Discussion of model results 

The model results for the outer lobe pair are broadly con- 
sistent with the properties of other old radio galaxies with 
an FRII-type morphology. The middle and inner lobes ap- 
pear to evolve inside the outer lobes and we describe here 
qualitatively how this may take place (see also Fig. 5). 

While the outer lobe is young and being inflated, it is 
overpressured with respect to the ambient medium. Thus the 
lobe drives a bow shock into the ambient medium, compress- 
ing and displacing it. Once the jet has stopped supplying 
energy to the lobe, the lobe gradually comes into pressure- 
equilibrium with the ambient medium. Those regions of the 
lobe in pressure equilibrium are no longer surrounded by a 
bow shock and they represent simply a bubble of light gas 
(lobe) immersed in the heavier ambient medium at the same 
pressure. The pressure balance is initially achieved close to 
the core where the pressure in the ambient medium is high- 
est; the central pinching of the lobe first discussed by Bald- 
win (1982). Note that the outer regions of the lobe may still 
be overpressured at this stage and therefore drive a bow 
shock. The parts of the lobe closest to the centre start to 
rise outwards due to buoyancy forces and are replaced by the 
heavier ambient medium sinking towards the core (labeled 
"Gas inflow" in Fig. 5). 

Buoyant replacement of the outer lobes of B0925+420 
sets in when the pressure inside these lobes becomes com- 
parable with the thermal pressure in the lobe environment. 
We can calculate a maximum pressure of 9 x 10 -15 J m -3 in- 
side the outer lobes allowed by the model parameters found 
in the previous sub-section. If we assume isothermal con- 
ditions in the lobe environment with a gas temperature of 
typically 10 7 K, then the thermal pressure of the ambient 
gas exceeds that inside the outer lobes at least out to a dis- 
tance of 25 ao ~ 50kpc. The buoyant replacement of the 
lobe proceeds at the sound speed of the ambient gas, about 
400 tons" for the assumed temperature. It then takes about 
40 Myr to buoyantly replace the outer lobe out to 20 kpc 
which is comparable with the size of the inner lobes. The 
same process would take roughly 200 Myr to replace the 
outer lobe out to 100 kpc. This is still considerably closer 
to the centre than the end of the middle lobes. Given the 
age of the outer lobes it is possible that the environment of 
the inner lobes is created by the buoyant replacement of the 
outer lobes. However, the same process is unlikely to be fast 
enough to provide the environment for the middle lobes. 

The middle lobes will initially have to plough their way 
through the ambient gas that has replaced the outer lobes. 
They may then reach the lower boundary of the outer lobes 
and continue to expand inside them. However, the modelling 
results suggest that (a) the middle lobes have also been re- 
placed by buoyant motions, and (b) the inner lobes are cur- 
rently expanding in that new ambient medium. The poor 
fit of our model to the observed properties of the middle 
lobes suggests that the radio synchrotron emission may not 
be that of the jet-inflated middle lobes themselves. 

5.5 Alternative model for the middle and inner 
lobe pairs 

Once the jets responsible for the middle lobe pair enter the 
outer lobes, they are immersed in a low density environment. 
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Clark & Burns (1991) simulated a restarting jet that enters 
the lobe it inflated in a previous activity phase. They find 
that the jet propagates almost ballistically through the old 
lobe and it does not form a significant new lobe inside the 
old one. However, they find a narrow bow shock driven by 
the jet into the material of the old lobe. This bow shock 
compresses the lobe material and may also re-accelerate rel- 
ativistic particles. The observed middle lobes of B0925+420 
are narrow and in the following we interpret them as the 
bow shocks surrounding the jets within the old lobes. 

We can derive the properties of the magnetised plasma 
in the outer lobes by assuming minimum energy conditions 
in these lobes (e.g. Longair 1984). For simplicity we assume 
that the power-law energy distribution of relativistic elec- 
trons, in terms of the electron Lorentz factor, extends from 
7 = 1 to 7 — > oo with an exponent of 2.5. We also assume 
a cylindrical geometry for the outer lobes with the emitting 
region at 326 MHz extending about three quarters of the 
way back from the leading edge towards the source centre. 
We then derive the magnetic field strength for both outer 
lobes as B; = 0.2 nT, where i refers to the initial, or pre- 
shocked plasma. The integrated density of electrons in the 
assumed energy distribution is then m = 0.09 m~ 3 . We use 
the same assumptions for the middle lobes, except that the 
radio emission of these structures observed at 1.4 GHz only 
extends to about halfway between the leading edge of each 
lobe and the source centre. We find a magnetic field strength 
of B a — 1 nT for both middle lobes (i.e. the shocked plasma) 
which implies an electron density of around n s — 2 m" 3 . 

If the radio emission of the middle lobes arises from the 
shock-compressed material in the outer lobes, then B s and 
n s must be related to £?; and m by the relativistic shock 
jump conditions (e.g. Landau & Lifshitz 1987). We can test 
this by deriving a value for the shock velocity (in units of c, 
the speed of light) in the rest-frame of the unshocked plasma 
in the outer lobe from 



(Ps - Pi) (e s +pi) 
(e s - e;) (a +p s ) ' 



(2) 



where Pi,p 3 are the plasma pressures measured in the 
respective rest-frames of the initial and shocked plasma re- 
spectively. The energy densities ei,e a are defined as e = 
nine 2 + u, with m the average rest-mass of a gas particle 
and u the thermal energy density, again measured in the 
rest-frame of both the initial and the shocked plasma. For 
simplicity we assume that the outer lobes contain either an 
ultra-relativistic pair plasma with u = 3p and m = m or a 
proton-electron plasma with one cold proton for each rela- 
tivistic electron with u — (3/2)p and m = m p . 

In both cases the pressure terms are dominated by the 
magnetic field and the relativistic particles and it then fol- 
lows from the assumption of minimum energy conditions 
u = 7B 2 /(6/io), with fio the permeability of the vacuum. 
For the case of a pure pair plasma we find the unphysical 
result of f3 B > 1. Clearly the mass density of a pure pair 
plasma is too low to be consistent with the shock jump con- 
ditions. In other words, a pure plasma in the outer lobes 
is insufficient to stop the jets and force the formation of a 
bow shock within the outer lobes. On the other hand, for 
the case of a proton-electron plasma we find (3 B = 0.3. The 
propagation of the bow shocks through the outer lobes is 



very fast. The jets are therefore almost ballistic as predicted 
by the numerical simulations. 

We can also directly study the momentum balance be- 
tween the jet and the receding material in the outer lobe. 
In the rest-frame of the bow shock, momentum balance re- 
quires (Bicknell 1994) 

2 a 2 



(3) 



where Wj and pj are the enthalpy and the pressure in 
the jet with a cross-section Aj. For a jet dominated by rel- 
ativistic particles Wj ~ 4pj. A B is the area of the working 
surface of the bow shock, while f3 B j is the velocity and 7 s j 
the Lorentz factor of the jet material measured in the rest- 
frame of the bow shock. It can be shown that the relativistic 
addition of velocities results in 



7sj/?sj = 7s 2 Pf (A - /3 S ) 2 : 



(4) 



where f3j and 7j are the velocity and Lorentz factor of the jet 
material in the rest-frame of the unshocked lobe material. 
Therefore 

[4% 2 7 2 (A - A) 2 + 1] PiA } = n t m(3 2 B c 2 A B . (5) 

The jet power of a jet dominated by relativistic par- 
ticles, i.e. the ratio of rest-mass energy density to internal 
energy of the jet material, 1Z, approaches zero, is (Bicknell 
1995) 



Q = 4cA iPn 2 f3 i . 



(6) 



Solving for p s A s and substituting into the equation for 
momentum balance we find 



_ 4 7s 2 7 2 (A-A) 2 + l 



(7) 



This expression for the area of the working surface of 
the bow shock is not very sensitive to the velocity of the jet 
material. For a jet power Q — 10 38 W, reasonable values of 
7j < 10 and again assuming a proton-electron plasma in the 
outer lobes, we get A s ~ 100 kpc 2 corresponding to a radius 
of R B ~ 6 kpc. This is in good agreement with the observed 
width of the middle lobes at their leading edges. 

We can apply the same arguments to the inner lobes and 
interpret them as bow shocks inside the outer lobe structure. 
In principle the inner lobes may in fact be located inside 
the middle lobes. However, if the middle lobes are indeed 
bow shocks, then they will not form wide lobes and hence 
it is more likely that the inner lobes are also located inside 
the outer, rather than the middle, lobes. In this case, the 
magnetic field strength inside the inner lobes is 2nT and 
the electron density is 9 m -3 . A pure pair plasma can again 
be ruled out, because it leads to a bow shock advance speed 
in excess of c. For a proton-electron plasma we find f3 B = 0.5 
and R B ~ 3 kpc. The inner lobes are consistent with both 
model interpretations. 

For the interpretation of the inner and the middle lobes 
as jet-driven bow shocks within the outer lobes we require 
a gas density in the outer lobes equivalent to at least one 
proton per relativistic electron. This is a lower limit and 
the real density in the outer lobes may be higher. Clearly 
the protons may have been transported into the outer lobes 
by the jets that inflated them. However, this would signif- 
icantly increase the already large power of these jets esti- 
mated in the previous sub-sections on the basis of the radio 
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luminosity of the outer lobes alone. The alternative is that 
the protons have "contaminated" the lobes across their sur- 
face. Kaiser et al. (2000) found that fluid instabilities on the 
lobe surface are slow to develop sufficiently to lead to signif- 
icant mixing across the lobe boundary. Instead, formation 
of an atmosphere of sufficient density for the ram pressure 
confinement of the younger jets could take place via the dis- 
persion of warm gas clouds. These clouds are passing from 
the ambient medium through the bow shock into the lobes 
where they are slowly disrupted by turbulent motions (e.g. 
Gregory et al. 1999). As suggested by Kaiser et al. (2000), 
the gas stripped off warm clouds inside the lobes may pro- 
vide a dense enough medium to confine the middle lobes of 
B0925+420. 



6 COMPARISON WITH GALACTIC X-RAY 
TRANSIENTS 

There is great interest, currently, in comparison of AGN and 
Galactic black hole candidates, in particular those Galactic 
transients which vary rapidly, typically showing outbursts 
which last a few months (e.g. Brocksopp, Bandyopadhyay 
& Fender 2004; Brocksopp et al. 2006). Correlations be- 
tween the X-ray and radio emission have been found, ex- 
tending over many orders of magnitude in both luminosity 
and mass (Merloni, Heinz & di Matteo 2003; Falcke, Kording 
& Markoff 2004). Furthermore there are suggestions that 
certain types of AGN are analogous to phases of jet/disc ac- 
tivity in the Galactic transients (Kording, Jester & Fender 
2006). For the Galactic transients, jet activity is very closely 
linked with the X-ray spectral behaviour of the accretion 
disc and corona mass. Repeated interactions with a com- 
panion galaxy was postulated as a mechanism by which the 
accretion rate could similarly be varied (Schoenmakers et 
al. 2000a). Of course, the timescales for variability in the 
AGN are just too long for us to make a direct comparison 
between AGN and Galactic transients. The DDRG provide 
a valuable link in this picture; while they still vary on unob- 
servable timescales, we are now able to obtain evidence for 
previous activity and thus obtain a more accurate idea of 
how the variability of jets in AGN and Galactic transients 
can be compared. 

A further similarity between the Galactic transients and 
some lobes of some DDRG is the absence of hotspots (e.g. in 
the outer lobes, and perhaps the middle lobes, of B0925+420 
- see Section 2). Very few of the Galactic objects have had 
extended emission from lobes or shocks detected, with the 
exception of e.g. Cyg X-l (Gallo et al. 2003), 1E1740. 7-2942 
(Mirabel et al. 1992) and perhaps GRS 1915+105 (Kaiser 
et al. 2004). Galactic sources and the DDRG have powerful 
radio jets, the energy of which needs to be dissipated in the 
ISM/IGM. If this is not through shocks at the head of the 
jet (i.e. hotspots) then it is not clear exactly how this dissi- 
pation takes place. It is particularly difficult to study this in 
the Galactic sources since it is rare that the jet emission is 
resolved; observations of both the jet and the bow shock in 
Cyg X-l have led to suggestions that the jet energy may be 
dissipated via the transport of protons (Heinz 2006). For the 
DDRG, it seems more likely that, over the large timescales 
involved, the jets have simply stopped providing energy to 



the lobes and any residual energy can be gradually dispersed 
in the IGM. 



7 CONCLUSIONS 

We present radio images of B0925+420, the first known 
"triple-double" FRII radio galaxy and resolve, for the first 
time, the inner pair of radio lobes. B0925+420 provides the 
strongest evidence yet that the class of DDRGs represent 
multiple episodes of jet activity and are not just co-aligned 
knots within a fainter jet. We also image this galaxy in lin- 
ear polarisation and find that the lobes are polarised up to 
~ 20%. We model the flux densities and sizes of the lobes 
to determine permitted ranges of jet power in terms of the 
age of the source and density of the surrounding medium. 
The outer and inner lobes are consistent with the model, 
and we conclude that the inner lobes lie within the cocoon 
of the original jet. The middle lobes are more complicated, 
requiring higher densities than those found within the outer 
lobes. An alternative model, interpreting the middle and in- 
ner "lobes" as additional bow shocks within the outer lobes, 
may be more successful. We suggest that some additional gas 
has refilled the region, perhaps via the dispersion of warm 
gas clouds. 
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